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ABSTRACT 

The stellar halos of spiral galaxies bear important chemo-dynamical signatures of 
galaxy formation. We present here the analysis of 89 semi-cosmological spiral galaxy 
simulations, spanning ^ 4 magnitudes in total galactic luminosity. These simulations 
sample a wide variety of merging histories and show significant dispersion in halo 
metallicity at a given total luminosity - more than a factor of ten in metallicity. Our 
preliminary analysis suggests that galaxies with a more extended merging history 
possess halos which have younger and more metal rich stellar populations than the 
stellar halos associated with galaxies with a more abbreviated assembly. A correlation 
between halo metallicity and its surface brightness has also been found, reflecting the 
correlation between halo metallicity and its stellar mass. Our simulations are compared 
with recent Hubble Space Telescope observations of resolved stellar halos in nearby 
spirals. 
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1 INTRODUCTION 

Understanding the formation history of stellar halos is one 
of the classical pursuits of galactic astronomy. The prob- 
lem is generally framed within the context of two com- 
peting scenarios: one of "rapid collapse" (Eggen, Lynden- 
Bell & Sandage 1962), in which the stellar halo is formed 
by the rapid collapse of a proto-galaxy within a dynam- 
ical timescale (~ 10 ^ yr), and one of "galactic assembly" 
JSearle fc ZinnI Il978l) . whereby the stellar halo is assem- 
bled on a longer timescale (~ 10^ yr) by the accretion of 
"building-blocks", each with separate enrichment histories. 
Both scenarios have their strengths and weaknesses, and it 
would appear that a hybrid model is t he most plausible op- 
tion consistent with extant data (e .g.: IChiba fc Beer jf200(i 
iFreeman fc Bland-Hawthornll200 . 

An intriguing piece of the halo formation "puzzle" is 
provided by comparing the stellar halo of our own Milky 
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Way with that of its neighbour, M31. First, despite their 
comparable total galactic luminosities, the stellar halo of 
M31 is apparently much more metal-rich than that of the 
Milky Way (e.g.: iR.van fc Norri3 Il99ll: iMould fc KristianI 
198d: Durrell, Harris fc Pritchet 2001: iFerguson et al.lj 2002r 
Brown et al.l20 0g':'Bellazzini et al.'2003l: lDurrell et alj2004l : 
Ferguson et ab 2005: Irwin ct al. 2005). In fact, the halo of 
M31 bears a closer re semblance to that of NGC 5128 (e.g., 
iHarris fc Harri!j|200J) . despite their differing morphological 
classifications. 

The stellar halo - galaxy formation symbiosis has been 
further brought to light by the recent work of Mouhcine et 
al. (see also Tikhonov, Galazutdinova fc Drozdovsky 2005). 
The deep Hubble Space Telescope (HST hereafter) imag- 
ing of nearby spiral galaxy stellar halos in lMouhcine et all 
[2005^ suggests a significant correlation between stellar halo 
metallicity and total galactic luminosity. On the surface, this 
correlation appears to leave our own Milky Way's halo as 
an outlier with respect to other spirals, with a stellar halo 
metallicity ~1 dex lower than spirals of comparable lumi- 
nosity (such as M31, as alluded to earlier). 
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However, it must be noted that the metallicity of the 
Galactic halo in the Solar Neighbourhood (e.g., in Ryan 
& Norris 1991) comes from spectroscopic metallicities in 
a kinematically-selected sample, whereas that of the M31 
halo, as well as those of the stel lar halos of nearby spi- 
ral galaxies (|Mouhcine et alJl2^^!^^ . have been derived pri- 
marily from photometric metallicities in geographically- 
selected samples. Because it is now becoming possible 
to obtain spectroscopic metallicities for significant sam- 
ples of kinematically-selected giants in the M31 halo (e.g.: 
Ilbata et al]l2004l: ICuhathakurta et alj|2005|) . it wiU be cru- 
cial to assess the consistency of the M31 stellar halo metallic- 
ity as derived from kinematically- and from geographically- 
selected samples, respectively. 

A grasp of the true scatter around the general trend 
in the halo metallicity - luminosity relation awaits a larger 
observational data set, however any theory which attempts 
to explain such a relation needs to simultaneously account 
for the apparent metallicity discrepancy between the Milky 
Way and the Andromeda stellar halos. The question arises as 
to what is driving the scatter of halo metallicity for galaxies 
of comparable luminosity? 

In what follows, we investigate whether difi'erences in 
the stellar halo assembly history can explain the diversity 
seen in halo metallicities. Using chemo-dynamical numeri- 
cal simulations, we have constructed a sample of 89 model 
disc galaxies, spanning ~4 magnitudes in luminosity, and 
sampling a wide range of assembly histories at a given lumi- 
nosity (or mass). We contrast the derived stellar halo metal- 
licity - galactic lum inosity relation with th e recent empiri- 
cal determination of iMouhcine et al] f2005) . The numerical 
framework in which the simulations have been conducted is 
described in Section 2, while Sections 3 and 4 present the 
results and the related discussion, respectively. 



2 SIMULATIONS 

The simulations employed here are pattern ed after the 
semi- cosmological adiabatic feedback model of lBrook et alJ 
i2004l) . an d were constructed using the chemo-dynamical 
code GCD+ iKawata fc Gibsonll2003l) GCD+ self-consistently 
treats the effects of gravity, gas dynamics, radiative cooling, 
star formation, and chemical enrichment (including Type II 
and la Supernovae, relaxing the instantaneous recycling ap- 
proximation) . 

This semi-cosmological version of GCD+ is based on the 
seminal work of Katz & Gumi (199 J). The initial condi- 
tion for a given model is an isolated sphere of dark mat- 
ter and gas. This "top-hat" overdensity has an amplitude 
5i at initial redshift Zi, which is approximately related to 
the collapse redshift Zc by Zc = 0.365i(l + Zi) — 1 (e.g., 
fPadman abhan 1993il . We set Zc = 2.0, which determines 
5i at Zi = 40. Small-scale density fluctuations based on a 
CDM power s pectrum are super imposed on the sphere us- 
ing COSMICS^ iBert8chingeilll998l'l . and the amplitude of the 
fluctuations is parameterised by (Tg. These fluctuations are 
the seeds for local collapse and subsequent star formation. 
Solid-body rotation corresponding to a spin parameter A is 
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imparted to the initial sphere to incorporate the effects of 
longer wavelength fluctuations. For the flat CDM model de- 
scribed here, the relevant parameters include Oq = 1, baryon 
fraction Q,b — 0.1, -ffo = 50 km s^^ Mpc~^, spin parameter 
A = 0.06, and as, — 0.5. We employed 14147 dark matter and 
14147 gas/star particles, for the 89 simulations presented 
here. 

Using the same number of particles, collapse redshift 
Zc, and spin parameter A, a series of 89 simulations were 
completed spanning a factor of 50 in mass in four separate 
mass "bins": Mtot = IxlO" Mq; 5xl0" Mq; 1x10^^ M©; 
5x 10^^ M0. For each total mass, we run models with differ- 
ent patterns of small-scale density fluctuations which lead 
to different hierarchical assembly histories. This was con- 
trolled by setting different random seeds for the Gaussian 
perturbation generator in COSMICS. 



3 RESULTS 

Our grid of 89 simulations were employed to populate the 
redshift z = Q stellar halo metallicity - luminosity plane 
(i.e. ([Fe/H])haio - My) shown in Fig. 1 (flUed symbols^). 
We have been conservative in our geometrical deflnition of 
"stellar halo", adopting a projected cut-off radius of 15 kpc 
to delineate the halo from (possible) disc-bulge contami- 
nants. Such a choice should also allow an easier comparison 
with the observations of halo flelds in nearby spiral galaxies. 

The stellar metallicity distribution function (MDF) was 
then generated using all stellar particles from the halo "re- 
gion", convolved with a Gaussian of cr[Fe/H] = 0.15 dex, 
representing the typ ical observational uncertainties (e.g., 
iBellazzini et al.ll2003ll . The peaks of the derived MDFs are 
shown^ in Fig. 1; 68% of the stellar particles within a given 
(simulated) MDF typically span 1.4 dex in metallicity. We 
confirmed that the derived halo MDFs are insensitive to 
projection effects, even in the most massive (largest) mod- 
els, suggesting that our i? ~ 15 kpc definition for the halo 
"region" is sufficient for minimising disc and bulge contam- 
ination. The optical properties of our simulated stellar pop- 
ul ations were derived us ing the population synthesis models 
of iMouhcine fc LanconI (2003), taking into account the age 
and the metallicity of each stellar particle. We note that the 
V-band luminosity shown along the abscissa of Fig. 1 is the 
total galactic luminosity (i.e. halo -I- bulge -I- disc). 

Fig. 1 also shows the corresponding observed halo 
metallicity - lu minosity values f or 13 nearby spirals, taken 
primarily from IMouhcine et a l. ('2005') , supplemente d with 
data from the literature iBinnev fc Merrifield. .1998 : lElsonI 
I1997I: Harris, Harris & Poole 1999; Brooks, Wilson & Har- 
ris 2004)''. For the observational datasets, we calculated the 
MDF peak and 68% "dispersions" exactly as we did for the 



^ Note: MDF centroids were constrained to sit on a 0.2 dex "grid" 
in ([Fe/H])iiaioi which is the source of the apparently "quantised" 
centroid values in Fig. 1. 

^ Only those simulations with >100 halo stellar particles are in- 
cluded in the analysis here. 

^ Note that the identification of the stellar halo in M33 is still 
debated (e.g., Tiede, Sarajedini & Barker 2004). 
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Figure 1. The stellar halo metallicity — total galactic V-band lu- 
minosity relation (([Fe/H])ijaio ^ My)- Filled squares correspond 
to the peak of the MDF for each of our simulations. Open (un- 
labelled) squares represent the data of Mouhcine et al. (2005); 
open (labelled) squares represent additional data taken from the 
literature (see the text for details). 68% of the stars in the sim- 
ulated (observed) MDFs are typically enclosed within ±0.7 dex 
(±0.35 dex) of the peak of the respective MDF. 

simulated datasets^; 68% of the stars in the observed MDFs 
are typicaUy enclosed within ±0.35 dex of the MDF peak, a 
factor of ~2 narrower than the simulations^ . 

What is readily apparent from Fig. 1 is that significant 
variation in halo metallicity (>1 dex) exists at any given to- 
tal galactic luminosity in our simulations. For a set of models 
with the same initial total mass, the only difference among 
these models can be traced to the random pattern of ini- 
tial small-scale density fiuctuations; this translates directly 
into differing hierarchical assembly histories. Qualitiatively, 
it would appear that assembly history alone may account 
for the diversity in halo metallicity at a given galactic lu- 
minosity, and thus account for the apparent outliers in the 
observed trend. 

Fig. 2 shows snapshots of the gas particles at redshift z 



^ Each MDF has been fitted by an u nivari- 
ate skew-normal distribution (e.g.: [Azzalin? '2 0051 
|http : //azzalini . Stat .unipd. it/SN/Intro/intro .html I. For 
the observational datasets, we calculated the MDF as derived 
from the mean location of the Red Giant Branch stars in 
the Colour Magnitud e Diagrams of the fields observed in 
iMouhcine et alj i2005l) . 

^ Our current model does not take into account any pre- 
enrichment scenario due to extremely metal— poor stars (Pop III 
hereafter) whose detailed phys ics is still m uch debated (e.g.: 
Woosley, Heger & Weaver 2002; |Lareo3|20o3). An early and ho- 
mogeneous pre— enrichment of the proto-galactic masses due to 
Pop III could lead to narrower MDFs at lower redshift. 



Figure 2. Upper panels: Projected distribution of the gas parti- 
cles at redshift z = 1.5 for the simulations with Mtot = 10^^ Mq. 
Lower panels: The associated stellar age distributions at z = 
for the galaxies in the upper panels. The solid (dotted) histogram 
corresponds to the stellar age distribution function for the entire 
galaxy (stellar halo). The corresponding halo metallicities are de- 
noted in each panel. 

= 1.5 and the associated z = stellar Age Distribution Func- 
tion (ADF) for three of the models with Mtot = 10^^ Mq. 
These representative models demonstrate that the simulated 
galaxies with the more metal-rich halos are assembled over 
a longer timescale, and thus possess a broader ADF (both 
in the halo and in the associated galaxy). Su ch a scenario 
is consistent with the evidence, presented by iBrown et alJ 
(2003), of a substantial intermediate-age metal-rich popu- 
lation in a geographically-selected halo field of M31^. Con- 
versely, the simulated galaxies with the more metal-poor 
halos are assembled earlier through more of a monolithic 
process, with a consequently narrower ADF. Our simula- 
tions suggest that the halo metallicity refiects directly the 
formation and assembly history of the host galaxy*. 

Further, the most metal-poor stellar halos in our sim- 
ulations (which, recall, formed preferentially via more of 
a monolithic collapse) possess a-elements to iron ratios a 
factor of ^^2 higher than the most metal-rich halos (which 
formed preferentially over more extended period of hierar- 

However, the eventual contamination (by the disc or b y debris 
of accrete d satellites) in th is outer field is debated (e.g.: ISrownl 
2003; Ferg uson et alil2005l) . 

We note that the "dispersion" in the halo and total galactic 
ADFs appear correlated, such that galaxies with narrower halo 
ADFs also appear to have narrower total galactic ADFs. This 
could suggest that the assembly history of the halo "knows" of 
the formation history of the other galactic structural components 
(e.g., bulge and disc). We will explore this suggestion in a future 
study. 
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Figure 3. The relation between the peak of the halo [0/Fe] distri- 
bution function and the peak of the halo MDF in our simulations. 



chical clustering), as shown in Fig. 3. Such a trend is ex- 
pected if the different amount of a-elements and Iron re- 
leased by Type II and la Supernovae over different timescales 
is taken into account (e.g.; Timmes, Woosley & Weaver 1995; 
IWooslev et al]|2002l) . 

The upper panels of Fig. 4 show the stellar halo metal- 
licity - I-band surface brightness relation (([Fe/H])haio - 
/ij"*'") for models with Mtot = 10^^ Mq; the surface bright- 
ness was measured at a projected distance of 20 kpc from 
the dynamical centre of each simulation. The halo metallic- 
ity - stellar mass relation ({[Fe/H])haio - M1^^'°) is shown in 
the lower panels of Fig. 4. The three galaxies presented in 
Fig. 2 are also labelled in Fig. 4. An immediate correlation 
is apparent with the more massive halos possessing higher 
surface brightness and also higher metallicity. This is consis- 
tent with a picture in which galaxies that experienced more 
extended assembly histories have more massive stellar halos, 
with both higher halo metallicities and halo surface bright- 
nesses - i.e. higher stellar halo densities. Similar trends are 
observed in our other total mass models, as shown in the 
right panels of Fig. 4. 

The issue of model convergence (resolution) is always a 
concern when interpreting cosmological simulations (partic- 
ularly when including baryons). To test this, we conducted 
a series of simulations for Mtot =10^^ Mq, with 5575x2 and 
9171x2 particles, to supplement the default grid (which used 
14147x2 particles); the middle panels of Fig. 4 show where 
these lower-resolution simulations sit in the halo metallicity 
- surface brightness plane and in the halo metallicity - stel- 
lar mass plane, respectively. The consistency between higher 
(left panels) and lower (middle panels) resolution models is 



Figure 4. Upper panels: Stellar halo metallicity - halo I-band 
surface brightness relation ({[Fe/H])iiaio - ^5"''°): in the left panel, 
our fiducial models with Mtot = 10^^ Mq and 14147x2 parti- 
cles are shown, while the middle panel shows the corresponding 
lower-resolution models with 5575x2 (crosses) and 9171x2 (open 
squares) particles. The metallicity error bars indicate the 68% 
Confidence Level. The average number of halo particles is shown 
in the bottom right corner of each panel. Models with fewer than 
100 stellar halo particles (R>15 kpc) are not included here. Mod- 
els A - C of Fig. 2 are denoted in the left and middle panels. 
The right panel shows the ([Fe/H])halo - ^t^'''" relation for all the 
simulations. Lower panels: Halo metallicity — stellar mass relation 
(([Fe/H])iiaio - Mj;^'°): in the left panel, our fiducial models with 
Mtot = lO^'^ Mq, while the middle panel shows the corresponding 
lower-resolution models. The right panel shows the ([Fe/H]}ijaio 
— M^^'° relation for all the simulations. 

reassuring and leads us to conclude that our simulations are 
not affected significantly by numerical resolution.® 



4 SUMMARY AND CONCLUSIONS 

We have presented here an analysis of the characteristics of 
spiral galaxy stellar halos formed within a large grid of nu- 
merical simulations, with particular emphasis placed upon 
the relationship between stellar halo metallicity and the as- 
sociated galactic luminosity. We have demonstrated that at 
any given total luminosity (or, conversely, total dynamical 
mass), the metallicities of these simulated stellar halos span 
a range in excess of 1 dex. We suggest that the under- 
lying driver of this metallicity spread can be traced to the 
diversity of galactic mass assembly histories inherent within 

^ Because we limit our analyses to models with >100 halo stellar 
particles, the /ij?'''°>28 mag arcsec"'^ region of the halo metal- 
licity - surface brightness plane and the MJJ'''°<10®Mq region of 
the halo metallicity - stellar mass plane are underpopulated by 
the lower— resolution models (middle panels of Fig. 4). 
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the hierarchical clustering paradigm. Galaxies with a more 
protracted assembly history possess more metal-rich and 
younger stellar halos, with an associated greater dispersion 
in age, than galaxies which experience more of a monolithic 
collapse. 

For a given total luminosity (or dynamical mass), those 
galaxies with more extended assembly histories also pos- 
sess more massive stellar halos, which in turn leads to a 
direct correlation between a stellar halo's metallicity and its 
surface bright ness (as anticipate d by earlier semi-analytical 
models - e.g.. lRenda et al]l2005ll . By extension, such a cor- 
relation may prove to be a useful diagnostic tool for disen- 
tangling the f ormation history of dis c galaxies. 

Recentlv. iMouhcine et all ^2005^ have presented an ob- 
served correlation between stellar halo metallicity and total 
galactic luminosity, as shown in Fig. 1. The observed disper- 
sion in the mean halo metallicity at a given galactic luminos- 
ity is smaller than what we find in our simulations. However 
the latter can account for the outliers in the observed trend. 
Since our motivation has been to study which is the effect 
of the pattern of the initial density fluctuations alone on 
the stellar halo features at redshift zero in simulated spiral 
galaxies, it is worth to note that galaxy formation, as it is 
observed, is an ongoing process which is the result of the 
interplay among different parameters, of which the pattern 
of initial density fluctuations (thus the merging history) is 
one. We have shown that the merging history alone may 
be held responsible of the dispersion in halo metallicity at a 
comparable total galactic luminosity, as apparently observed 
for example in our Milky Way and in Andromeda (see Sec- 
tion 1). 

This begs the question... FPTizc/i is normal? Our study 
suggests that if the stellar halo was Eissembled (primarily) 
through more of a monolithic collapse, such a low metallicity 
is indeed what should be expected; conversely, the fact that 
the M31 stellar halo is significantly metal-rich is suggestive 
of a more protracted assembly history. An observational con- 
sequence of these differing formation histories is the predic- 
tion that the M31 stellar halo should possess a high surface 
brightness; observations tentatively suppo rt this prediction 
(Reitzel, Guhathakurta & Gould 1998: lGuhathakurta et aD 
l2005l:llrwin et alJl2005^ . 

Further observations are needed to tighten our grasp 
of the strength and scatter of the stellar halo metallicity 
- luminosity relation, which we have shown to be a useful 
diagnostic tool for disentangling the formation history of 
disc galaxies. 
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